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Abstract- In this work, one uses the method of 
Monte Carlo to simulate a luminescent discharge in 
continuous, in a gas electropositive  (argon) a 
pressure of Itorr, this technique capable to couple 
the loaded particle movement (ion, electron) with 
the Poisson equation, it permit to determine the 
probability density function of the electronic 
energies (PDF) in plasma, the results are compared 
by its obtained by the kinetic and fluid models, [1], 
[2], they are similar. 
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Résumé — Dans ce travail, on utilise la méthode 
de Monte Carlo pour simuler une décharge 
luminescente en continu, dans un gaz électropositif 
(l’argon ) à une pression de Itorr, cette technique 
capable de coupler le mouvement des particules 
chargées (ion, électron) avec l’équation de Poisson, 
elle permet de déterminer la fonction de densité de 
probabilité des énergies électroniques (FDP) dans le 
plasma, les résultats sont comparés par les résultats 
obtenus par les différents modèles cinétiques et 
fluides ,[1],[2], ils sont similaires. 

Mots Clés : décharge luminescente, méthodes de 
Monte Carlo, décharge à basse pression. 


Ï) Introduction : 

L’axe des travaux dans les décharges 
luminescentes s’est tourné vers une modélisation 
plus détaillée de la décharge, une bonne 
modélisation permet de comprendre et de prédire 
leur comportement général, alors de contrôler et 
d'utiliser cette technologie dans les besoins 
industriels, telque, les semi-conducteurs, les 
panneaux de visualisation à plasma, le dépôt des 
couches minces, la gravure en micro-électronique et 
le traitement de surfaces[3]. Il existe trois types de 


modèles, le modèle fluide qui consiste à résoudre 
les équations hydrodynamiques[4], [5], [6], ces 
dernières utilisent les moments de la fonction de 
distribution, qui sont accessibles à la mesure, mais 
il n’est plus valable dans le cas de l’hors équilibre 
thermodynamique, le modèle cinétique basé sur la 
résolution directe de l’équation de Boltzmann [7], 
[8], ou bien par un autre modèle statistique basé 
sur des lois stochastiques et simule les trajectoires 
des particules chargées, la technique appelée 
méthode de Monte Carlo[9],[10]. 

ID) Description du modèle : 

Dans ce travail on utilise une géométrie 
cylindrique avec électrodes planes, le modèle à 
trois dimensions est réduit à deux dimensions 
axiale et radiale, dans ce modèle on considère 
qu’il n’y a pas d’effets radiaux (diffusion vers les 
électrodes), on considère aussi que le champ 
électrique est axial, et le modèle se réduit donc à 
une seule dimension, on garde la direction axiale 
et radiale de la vitesse. Le gaz de décharge est 
supposé à la température ambiante 300°k, on 
considère qu’il y’a quatre espèces des particules, 
les atomes neutres (Ar), les atomes excités (Ar'), 
les ions positifs (Ar*) et les électrons (e’), la 
densité des particules chargées est très inférieure 
devant la densité des neutres. Pour cette raison on 
s'intéresse seulement aux collisions entre les 
particules chargées et les neutres, ces collisions 
sont : 

-Les collisions élastiques, d’excitation et 
d’ionisation entre les électrons et les atomes 
neutres, 

-Les collisions de transfert de charge et les 
collisions élastiques entre les ions et les atomes 
neutres de l’argon, les atomes excités se 
désexcitent spontanément. 
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IT) Simulation des particules chargées par la 
méthode de Monte Carlo : 

La méthode de Monte Carlo consiste à simuler un 
ensemble d'évènements (particules chargées) par un 
autre ensemble facilement réalisable (nombres 
aléatoires). Elle consiste à considérer une 
population des particules chargées et retracer 
l’histoire de chacune d’elle depuis son émission 
jusqu’à sa disparition, les trajectoires des particules 
chargées sont données par la loi de Newton. 


Z=Zo +01 +5 E.Ai 4 


=" +V,oAt (1) 
Vz =v0+ LEA 


Vr =V, 
tel que zo et ro sont respectivement les positions 
axiale et radiale de la particule , z et r sont les 
positions axiale et radiale après un pas de temps A ; 
vo et vo Sont respectivement les vitesses axiale et 
radiale à l’instant #. v. et v, les vitesses axiale et 
radiale à l’instant fo+Af, E:0 est le champ électrique 
dans la position zo, et g et m sont respectivement la 
charge et la masse de la particule, l’énergie 
cinétique est donnée par la relation suivante : 


env, +, (2) 


2 
Le champ électrique est obtenu par la résolution 
de l’équation de Poisson : 


Te = Line) @ 
__ dU() 
Ne ra 


ni(z) et ne(z) sont respectivement la densité des ions 
positifs et la densité des électrons. 

La probabilité de collision, la nouvelle énergie et la 
direction juste après la collision sont déterminées 
par un tirage des nombres aléatoires. 
P=I-exp(-010(E)AS) (4) 
où P : la probabilité totale de collision, As est la 
distance entre le point (z0,ro) et le point (z,r), no la 
densité des neutres et Ge) la section efficace 
totale de collision en fonction de l’énergie & .Le 
tirage d’un nombre aléatoire permet de dire si 1l y’a 
une collision ou non, si R<P il y’a une collision, si 
R>P pas de collision, maintenant si R<P, il faut 
ürer un autre nombre aléatoire (R ) pour déterminer 
le type de cette collision, en subdivisant l’intervalle 
des nombres aléatoires en sous intervalles chacun 
correspond à une fraction de probabilité, elles sont 
calculées à partir des relations suivants[11] : 


Pour les électrons : 
Oro €) =Ooni(E)+Oexc(E) +O etat €) 


P _ Oioni €) 
ionis— 
© eto( €) 
Oexc(E 
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Pour les ions : 


Otoi €) —O'ran-ch ne E) FO etat €) 
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OÙ Pionis Pexe et Peelas Sont respectivement les 
fraction de probabilité d’ionisation, d’excitation et 
des collisions élastiques avec les atomes neutres, 
Piran-charge  €t  Pieus Sont respectivement les 
fractions de probabilité de transfert de charge et 
des collisions élastiques avec les atomes neutres. 
L'énergie de chaque particule et sa direction sont 
déterminer aussi par tirage des nombres 
aléatoires, après la collision la particule fait un vol 
libre et la procédure se répète après un pas de 
temps. 

III-1) Simulation des électron : 

Les électrons secondaires émis par la cathode 
sous l’effet de bombardement ionique ont une 
énergie distribue entre O0 et 6 ev avec un 
maximum de 4 ev et suivant un angle axial @ et 
un angle radial @ , [11], tel que : 
cosô,=R, 
po =2TR2 (7) 
R1 et R2 des nombres aléatoires tiré au hasard, la 
nouvelle énergie et la direction dépendent du type 
de la collision. 

a) Ionisation: dans le cas d’ionisation l’énergie 
total avant la collision est divisé entre l’électron 
primaire et l’électron secondaire. 
£0 — Eïonis 

2 
E pri E0 — EÉionis — Esec (8) 


Esec=Rjonis 


OÙ &ec, épri Et &o Sont respectivement l’énergie de 
l’électron secondaire, l’énergie de l’électron 
primaire et l’énergie de l’électron primaire juste 
avant la collision. Rionis est un nombre aléatoire et 
Eïonis l'énergie d’ionisation d’argon (&onis=15.76 
ev). 
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b) Excitation: la nouvelle énergie est donnée par la 
relation suivante : 

EE) —Eex (9) 
£exe et & sont respectivement l’énergie d’excitation 
(&xe—11.5 ev) et l’énergie de l’électron incident 
juste avant la collision. 

c) Collisions élastiques : l’énergie cinétique de 
l’électron après une collision élastique est donnée 
par la relation suivante : 
8=a [1-22 (1-05 D) 


at 


(10) 


OÙ Me €t Ma Sont respectivement la masse de 
l’électron et la masse de l’atome cible et 7 l’angle 


CA E 


de diffusion axial dans le repère de collision. 


fig.1 :schéma représentant la séquence des collisions. 


L’angle axial juste après la collision est donnée 
par la relation suivante : 
cos@=cos6, cos 7—sin 0, sin 7cosy (11) 
telque @ l’angle axial juste avant la collision, x et y 
les angles de diffusion axial et azimutal. Ces angles 
sont calculer à partir de tirage des nombres 
aléatoires R, et Ry. 
cos 7=I-R; 


y=27R (12) 


14 


10 


* collisions élastiques 
x l’ionisation 
o l’excitation 


fig.2 : les sections efficaces de collision e/Ar [12],[13] 


II-2) Simulation des ions positifs : 

Dans le cas d’ionisation non seulement un 
électron est produit, mais aussi un ion positif, le 
transport de cet ion est simulé par la méthode de 
Monte Carlo, l’indépendance de simulation entre 
les électrons et les ions est due à la faible 
ionisation du gaz (les collisions électron-ion sont 
négligeables). L’énergie et la direction de l’ion 
juste après la collision dépend avec le type de 
collision. 

a) Transfert de charge: dans ce type de 
collision la charge de l’ion se transfert vers 
l’atome et le nouveau ion démarre avec une 
vitesse nulle suivant la direction du champ, la 
section efficace de ce type de collision est 
supposée constante Oiran-charge =10°!* cm?{[14 ], 

b) Collisions élastiques : les collisions entre les 
ions et les atomes d’argon sont considérées 
comme des collisions entre des sphères dures, la 
section efficace est donnée par la relation 
suivante : 

Oi,elas =47RA? (13) 
où RA- est le rayon de l’atome de l’argon, l’énergie 
de l’ion juste après la collision est donnée par la 
relation suivante : 


Eion=E0 ion ICOS? 7) (14) 


OÙ &ion l'énergie de l’ion juste avant la collision 
et y l’angle de diffusion axial dans le repère de 
collision. 

IV) Résultats et discussion : 

A titre d'illustration sur l’utilisation de la 
méthode de Monte Carlo, on présente une 
application à une décharge continue luminescente 
pour un gaz d’argon à 1torr de pression, un espace 
inter-électrode de 2 cm et une tension continue sur 
l’électrode de gauche de (—-300 v) et l’autre est 
mise à la terre, le coefficient d’émission 
électronique secondaire (y= 0.12). Les résultats 
numériques sont obtenus avec un pas en temps 
At = 2.10 llsec et un pas spatial Az=0.04 cm. 

La figure (3) représente la densité des particules 
chargées (électrons et ion positifs), cette densité 
est très importante dans le plasma que dans les 
gaines car les électrons émis par la cathode sous 
l’effet de bombardement ionique gagnent 
l’énergie dans la gaine cathodique et excite le 
milieu d’une façon intense dans la lisière de la 
gaine. Le plasma est caractérisé par une neutralité 
électrique d’ou un champ électrique presque nul, 
par contre les gaines sont caractérisées par des 
charges d’espace positives qui sont responsables à 
la chute cathodique et anodique. 
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+ Densité électronique 
* Densité ionique 


la densité en cm-3 


la distance en cm 


fig.3 : la densité des particules chargées dans une décharge 
luminescente en continu 
Les figures (4) et (5) représentent respectivement, 
la variation spatial de champ électrique et la 
variation spatial de potentiel, la chute dans la gaine 
cathodique due à la densité des charges d’espace 
positives, dans le plasma, le champ électrique est 
pratiquement nul et le potentiel est constant, car la 
densité des électrons est égale à celle des ions 
positifs, dans la gaine anodique la chute est due à la 
grande différence de mobilité entre les électrons et 
les ions et donc le potentiel du plasma est toujours 
supérieur à celle dans les électrodes. 


le champ électrique en (V/cm) 


la distance en cm 


fig.4 : le champ électrique dans une décharge luminescente 
en continu 


le potentiel en (V) 


la distance en cm 


fig.s : le potentiel dans une décharge luminescente en 
continu 


Les figures (6) et (7) représentent respectivement 
la variation spatiale des densités des courants 
électronique et ionique et la variation spatiale de 
la densité de courant total, ce dernier est égal à la 
somme de deux courants ionique et électronique, 
dans la gaine anodique la densité des électrons est 
trop faible, alors le courant total est un courant 
ionique, dans le plasma la densité des ions est 
importante, mais leur courant est faible à cause de 
leur faible mobilité et le courant total est un 
courant électronique, par contre dans la gaine 
anodique le courant électronique est supérieur à 
celle dans le plasma, ce dernier compense le 
courant ionique qui du à la chute anodique, le 
courant total est pratiquement constant à tout 
endroit de la décharge avec une valeur de l’ordre 
de 0.0171 mA/cm°. 


[ densité de courant ionique 
* densité de courant électronique 


la densité de courant en mA/cm2 


la distance en cm 
fig.6 : le densité des courants électronique et ionique dans 
une décharge luminescente en continu 


la densité de courant total en mA/cm2 


la distance en cm 


fig.7 : la densité de courant total dans une décharge 
luminescente en continu 


La figure (8) représente la fonction de densité 
de probabilité des énergies électroniques dans le 
plasma(FDP), elle est composée de deux parties 
avec une énergie moyenne de 6.29 eV , partie 
corps pour les faibles énergies où se trouve le 
grand nombre des électrons, ces derniers 
gouvernent pour l’essentiel les phénomènes de 
transports et la partie queue pour les énergies 
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supérieures au premier seuil d’excitation, ils sont 
responsables des collisions inélastiques. 


FDP 


l'energie en (eV) 
fig.6 : la fonction de densité de probabilité des énergies 
électroniques dans le plasma (FDP) 
V) Conclusion : 

L’élaboration d’un modèle statistique basé sur des 
lois stochastiques permet de simuler les trajectoires 
des particules chargées par une méthode de type 
Monte Carlo, le couplage avec l’équation de 
Poisson donne un modèle auto-cohérent. De ce 
dernier, on peut prédire les caractéristiques 
générales de ce type de décharges et de détermine 
d’une manière indirecte la fonction de densité de 
probabilité des énergies électroniques, d’où la 
résolution de l’équation de Boltzmann, on peut 
aussi introduit facilement le second coefficient de 
Towsend et étudier le comportement des électrons 
et des ions dans les situations hors équilibre 
thermodynamique (gaines), alors la prédiction et le 
contrôle du flux ionique et le potentiel surtout au 
voisinage de la cathode où les ions positifs sont 
accélérés et bombardent le substrat avec des 
énergies qui dépendent à la chute de potentiel dans 
la gaine, la distribution en énergie des ions 
détermine le dégât d’anisotropie de gravure. 
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